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As additional neutron scattering experiments are performed on a variety of high temperature su- 
perconducting compounds it appears that magnetic incommensuration is a phenomenon common 
to all of the samples studied. The newest experimental results indicate that incommensurate peaks 
exist at momentum q = 7r(l, 1 ± 5) and tt(1 ± S, 1) in La2- x Sr x Cu04 (LSCO) and YBa2Cu307_5 
(YBCO). The dependence of 8 with hole doping appears to be similar in both materials. In addi- 
tion, new ARPES data for LSCO as a function of doping show that its Fermi surface is qualitatively 
similar to the one of YBCO, contrary to what was previously believed. Early theoretical attempts to 
explain LSCO and YBCO behavior usually relied on one- or three-band Hubbard models or the t-J 
model with electron hopping beyond nearest-neighbor and with different parameter values for each 
material. In this paper it is shown that using a one band Hubbard U-t-t' model with a unique set of 
parameters, U/t=6 and t'/t=-0.25, good agreement is obtained between computational calculations 
and neutron scattering and ARPES experiments for LSCO and YBCO. It is also shown that using a 
more negative t'/t will induce short-range magnetic incommensuration along the diagonal direction 
in the Brillouin zone, in qualitative disagreement with the experimental results. At the finite tem- 
peratures of the present Monte Carlo simulation it is also observed that in this model the tendency 
to incommensuration appears to be more related to the shape of the two-dimensional Fermi surface 
and the strength of the interaction, rather than to charge order. 



PACS numbers: 71.10.Fd, 71.18.+y, 74.25.Ha, 75.40.Mg 
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I. INTRODUCTION 

Neutron scattering experiments continue providing ex- 
citing information about the behavior of spin and charge 
degrees of freedom in the high T c cuprates. Studies on 
La2_ x Sr x Cu04 (LSCO) 0] have shown the existence of 
incommensurability in the spin channel near the com- 
mensurate position (tt,tt). In this case, the incommensu- 
rate peaks were found at Qs — 7r(l ± S, 1) and 7r(l, 1 ± S). 
The value of S increases linearly with doping in the range 
0.05 < x < 0.14, while for x > 6.14, S plateaus at 0.25. @] 
The intensity at Q 7 = 7r(l ± 5/2, 1 ± 5/2) was observed 
to be 0.18 times that at Qs for x=0.14 while the inten- 
sity at Q T = (-7T, 7r) was negligible compared with the 
background. [^| For some time neutron scattering experi- 
ments in YBa2Cu307_5 (YBCO) were less clear. Rossat- 
Mignod and coworkers detected magnetic fluctuations 
only at the commensurate position Q , while Tranquada 
et.al Q noticed possible incommensurate fluctuations. 
Measurements by Dai et al. [^) using a new position- 
sensitive detector bank indicate that incommensurability 
in the spin channel indeed is present in YBCO. Although 
originally incommensurability was detected only along 
the diagonal direction in the Brillouin zone, the latest re- 
sults indicate that it occurs at Qs as in the case of LSCO, 
and the dependence of 5 with doping appears to be sim- 
ilar in both materials. In addition, Bi2Sr2CaCu20s 
has been recently studied by Mook and Chakoumakos 
H and an incommensurate fluctuation that occurs be- 
low T r was found. This incommensuration was identified 



with a dynamic charge density wave because its scatter- 
ing intensity appeared to increase with increasing mo- 
mentum transfer. Scattering that could be described as 
magnetic has not been observed yet in this material be- 
cause the experimental technique used does not allow 
to reach high enough values of the momentum trans- 
fer. However, magnetic excitations should be present if 
a dynamic strip phase gives origin to the charge peaks. 
Evidence of incommensurability in the charge channel 
was also observed in Lai.6- x Ndo.4Sr x Cu04 || which led 
to speculations about a similar behavior in LSCO. This 
would indicate that the incommensuration in the spin 
channel may be due to the existence of charge-stripe or- 
der rather than to some kind of charge uniform spiral 
spin state or 2D Fermi surface effects. |l(]] 

While some of these experiments suggest that mag- 
netic incommensurability in the cuprates maybe due to 
charge-stripe order in which the orientation of the stripes 
is not material dependent, the theoretical understanding 
of this phenomenon is less clear. Short-range magnetic 
incommensurate correlations in the spin channel were de- 
tected early on in the Hubbard jll],[l2) and t-J |l3| mod- 
els. The split of the commensurate peak was observed 
to be qualitatively similar to the behavior in LSCO and 
YBCO. However, the dependence of 5 with doping did 
not reproduce the experimental data and, as it will be 
shown later, the experimentally observed relative inten- 
sities of the peaks at different points in the Brillouin zone 
are not well reproduced either. No particular order was 
reported in the charge channel in these models through 
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numerical analysis. 0,^5) However, since the proposal 
of dynamical microphase separation in the CuC>2 planes 
[PH , a reanalysis of these conclusions is needed. 

At the same time that the recent neutron scattering 
experiments were discovering new similarities between 
the magnetic properties of the different cuprates, also 
recent angular-resolved photoemission (ARPES) exper- 
iments performed on LSCO mapped its Fermi surface 
(FS) ]l7j at different values of doping unveiling interest- 
ing similarities among the qualitative FS shape of several 
high T c materials. 

Motivated by the new experimental results the aim of 
this paper is to revisit the U-t-t' model, exploring numer- 
ically the spin and charge channels to determine whether 
it is possible to obtain agreement with the new data for 
LSCO and YBCO using a unique set of parameters in 
the model. It will be concluded that this is indeed possi- 
ble. The paper is organized as follows: in Section II the 
U — t — t' Hamiltonian and the notation used are intro- 
duced. Results on magnetic and charge correlations, as 
well as on the shape of the Fermi surface, are presented in 
Section III while Section IV is devoted to the conclusions. 



II. THE MODEL 



<5n r = J2 a cj CT c riCr — (n). Here (n) = 1 — x is the average 
density of electrons. The brackets in Eqs. (2) and (3) 
refer to thermal averaging in the grand canonical ensem- 
ble which will be performed using the standard Quantum 
Monte Carlo (QMC) determinantal method. 

Before presenting our results let us discuss the be- 
havior of the spin and charge correlations in the non- 
interacting system (U/t=0), and also in the standard 
Hubbard model with t'=0. In the non-interacting sys- 
tem the spin and charge correlations are related through 
S(q) = iN(q). N(q) increases from zero reaching the 
value (n) at q = 2kp and remaining constant afterwards. 
On the other hand, in the Hubbard model it was ob- 
served that at low density, i.e., (n) < 0.5, S(q) peaks 
at q = 2kF while N(q) is suppressed at these momenta 
compared with the non-interacting system and it only 
peaks at q = (ir,n). O] This peak is due to the short 
range effective repulsion between particles. The behavior 
at higher densities is very different in the spin channel. 
At half-filling a sharp peak develops at and, with a 
substantially reduced intensity, it moves to Qs with dop- 
ing. While the intensity is maximum at Qs as in the 
experiments, the weight at is always larger than the 
one at Q 7 1 1 which is qualitatively incorrect compared 
with recent experiments as it will be shown below. 



The U-t-t' one band Hubbard model Hamiltonian is 
jiven by 



III. RESULTS 



H = - t z2 ( c L c j^ + h - c -)-f v ( c [ Cn)(T +h.c.) 



<in>,cr 



U]T(n iT - 1/2) (n u - 1/2) + 



(1) 



where Cj CT creates an electron at site i with spin projec- 
tion a, m CT is the number operator, the sum (ij) runs over 
pairs of nearest neighbor lattice sites, and the sum (in) 
runs over pairs of lattice sites along the plaquette diago- 
nals. U is the on site Coulombic repulsion, t the nearest 
neighbor hopping amplitude, t' the diagonal hopping am- 
plitude, and fi the chemical potential. In this work t will 
be set equal to 1. 

The static charge and magnetic structure factors N(q) 
and S(q) are defined by the relations 



N(q) 



(cm (5n r ), 



S(q)=]Te^ r <S^ 



(2) 



(3) 



Due to the well known "sign problem" it is very diffi- 
cult to perform Monte Carlo numerical studies at small 
hole doping, low temperatures and values of U/t > 4. 
This problem is exacerbated as the absolute value of t' in- 
creases. For this reason the numerical efforts will be con- 
centrated here on the study of the fixed density (n) = 0.7, 
(i.e., x=0.3) since for this (n) a good degree of control of 
the numerical results can be achieved. In addition, new 
experiments have been performed in LSCO at precisely 
x=0.3 providing information about incommensuration |2j] 
and the shape of the FS (Tt]]. Here results for U/t=6 on 
8x8 lattices will be presented. Due to the sign problem 
at this relatively large value of U/t, and using a finite 
diagonal hopping t', the temperature had to be fixed at 
T=0.25t. It is to be expected that the intensity of the 
observed features will increase at lower temperatures and 
the results presented here are, thus, lower bounds to the 
actual values. 

Note that, as part of our study, runs for other values 
of U/t were performed and qualitative differences with 
the results for U/t=6 were not observed. In particular, 
for t'/t = we found that for values of U/t as high as 
10 the magnetic incommensuration always occurs at Q$ 
rather than along the diagonal as predicted by mean- field 
calculations in the strong coupling regime. |18|| 



where {SnoSn r ) and (SqS*) are equal-time density- and 
spin-correlation functions, S^ = \ J2 a /3 c r,a°a /3 c r,/3j an d 
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A. Magnetic Incommensurability 

As a first step, the static structure factor will be cal- 
culated for several values of t'/t and comparisons with 
the experimental neutron scattering results will be made 
along several directions in momentum space. 

In Fig.l-a the peaks in 5(q), indicative of short range 
spin incommensurate tendencies in the U — t — t' model, 
are presented along the (0, 7r) — (2tt, tt) direction for val- 
ues of t'/t ranging from to -0.5. The figure shows that 
for all values of t'/t analyzed here, the peak in the struc- 
ture factor occurs at q = (37r/4, tt) and (57r/4,7r) which 
correspond to S — 0.25. This is in agreement with the 
experimental value for x=0.3 in LSCO. [|| However, a 
spline fit through the available data points (dashed line) 
suggests that the actual peak at |t'/t| > 0.3 occurs at a 
slightly larger value of S. Actually, when the structure 
factor is scanned along the diagonal direction as shown 
in Fig.l-b it is clear that the results with |t'/t| = 0.3 or 
larger do not fit the experimental data because the max- 
ima in 5(q) at q = (1 ± <5)tt(1, 1) = (3tt/4, 3tt/4) and 
(57r/4, 57t/4) have intensities which are approximately 
equal or higher than those at Qs [ [l9| , as it can be de- 
duced by comparing Fig.l-a with Fig.l-b. This is in dis- 
agreement with the experiments that indicate that the 
intensity at q = (1 ± S)n(l,l) should be indistinguish- 
able from the background. 

An important result is that at V —0 S(q) along the di- 
agonal direction, q x = q y , has a maximum at q = Q n 
as it can be observed in Fig.l-b (top curve). This be- 
havior is in disagreement with the experimental data for 
LSCO presented in Fig. 3 (closed circles) of Ref. ||, and 
for YBCO in Fig.l-c of Ref. J(| where a minimum is ob- 
served at along the diagonal. The experimental re- 
sults were obtained at x=0.15 for optimally doped LSCO 
and x=0.1 (i.e., S = 0.4) for YBCO showing that the 
qualitative behavior does not depend strongly on dop- 
ing. This indicates that the standard Hubbard model 
(t' = 0) does not describe the qualitative behavior ob- 
served with neutron scattering when S = 0.25. However, 
working with t'=-0.2 the fit along the diagonal in mo- 
mentum space (Fig.l-b) indicates that S(Q 7V ) now has 
become a local minimum, in qualitative agreement with 
the experiments. Then, the constraints on the relative 
intensity of the peaks provided by the experiments leave 
a small window of possible values of t'/t. If |t'/t| > 0.3 
the incommensurate peaks will appear along the diago- 
nal rather than at Qs (as discussed in the previous para- 
graph) and if |t'/t| < 0.2 the structure factor at Q^ would 
be a local maximum rather than a minimum along the 
diagonal direction. 

Additional experimental measurements || have been 
reported along the direction Qs[= 7r(l — 5,1)] — » Q 7 [= 
tt(1 - 6/2, 1 + 6/2)} -> Q s [= tt(1, 1 + S)]. Note that for 
5 = 0.25 the point Q 1 is given by (77r/8, 9ir/8) which 
is computationally accessible only on 16 x 16 or larger 
lattices if periodic boundary conditions are being used. 



Since such a large lattice size is beyond the capabilities 
of the QMC at finite hole density we obtained the values 
of 5(q) at the equivalent Q 1 points located at 77r/8(l, 1) 
and 97r/8(l, 1) with the help of a spline fit of the available 
data in Fig.l-b (dashed line). 
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FIG. 1. The static structure factor S(q) for U/t=6, 
(n) = 0.7 on an 8 x 8 lattice for T=0.25t and values of t'/t 
ranging from to -0.5. a) Along the (n, 0) — (tt, 2n) direction; 
b) along the diagonal direction; c) along the q y — Sir + q x 
direction. The dashed line indicates a spline fit. 

In Fig.l-c the structure factor is shown along the line 
q y = St: + q x which corresponds to the direction along 
which the experimental neutron scattering data shown 
by open circles in Fig. 3 of Ref. || were taken. In the 
figure it can be seen that the relative intensity of the 
numerical data at q x — 7r(l — 8/2) — 77r/8, (i.e., Q 7 ) 
and q x = 7r(l ± S) = 0.757T or 1.25tt, (i.e., Qs) is a func- 
tion of t'. Considering, as in the experiments, that the 
intensity for q = 7r(l ± S)(l, 1) (this point corresponds 
to q x = 37r/4 and 57r/4 in Fig.l-b) has to correspond to 
the background, then the relationship Q-y/Qs = 0.18 is 
satisfied for t'/t=-0.25. This is used as a guide because 
this relationship may be doping dependent but since the 
window in t'/t is so narrow, at most an error of the order 
of 0.05 in the estimation of t'/t is being made. 

The above analysis shows that a comparison of the 
present numerical data when the measured incommen- 
surability 6 is 0.25 agrees with the data for LSCO with 
the same S using U/t=6 and t'/t=-0.25. Notice that the 
position and the relative intensity of the peaks does not 
change too much with temperature according to the ex- 
periment. 
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B. Fermi Surface 

The next issue that will be addressed is whether the 
ratio t'/t = —0.25, fixed by the spin structure factor anal- 
ysis in the previous subsection, will fit other experimental 
data such as, for example, the shape of the Fermi surface 
recently obtained using ARPES for LSCO. @ The pos- 
sible shape of the FS will be determined by analyzing 
the momentum distribution n(q) which is calculated by 
Fourier transforming the one-electron Green's function, 

gu = W 

cr 

that is evaluated using QMC. 

The criteria used here to obtain the most probable lo- 
cus of the FS from numerical n(q) data are two: (a) find 
the values of q where n(q) w 0.5 Q and (b) find the val- 
ues of q where n(q) changes the most rapidly p0| . For the 
case of t'/t=-0.25 we have observed that both methods 
provide similar results in the regions close to the diagonal 
direction in the Brillouin zone, but substantial differences 
were observed close to the (0, tt) and (w, 0) points. While 
criterion (a) indicated a FS closed around (0,0), criterion 
(b) indicated a FS closed around Since criterion 

(b) provided similar results for the non-interacting case 
with (n) = 0.7 and t'/t=-0.25, a situation where it is 
known that the FS actually closes around (0,0), it was 
decided that criterion (a) would be more effective in this 
context. 

The FS obtained with procedure (a) is shown in Fig. 2- 
a. The closed circles indicate where n(q) w 0.5 and the 
continuous line is a 6th order polynomial fit of the points. 
The obtained FS is very similar to the non-interacting 
one (denoted by a dashed line in the same figure). It 
satisfies Luttinger's theorem within error bars and, ac- 
tually it is in excellent agreement with the experimental 
data for x=0.3 |17| shown with open squares in Fig. 2- 
a. Thus, using tyt = —0.25 good agreement has been 
obtained between the numerical results and two inde- 
pendent experiments (neutron scattering and ARPES) 
performed in overdoped LSCO. Since the obtained Fermi 
surface resembles closely the non-interacting one, the FS 
in the U/t=0 limit was calculated for x=0.1 and its shape 
compared directly with the experimental results at this 
density. As it can be observed in Fig.2-b the agreement 
is once again very good. According to Ref. @ the FS 
for LSCO at optimal doping is still centered about (tt, tt) 
and it is qualitatively similar to the one obtained for 
x=0.1. This is indeed what happens with the U/t=0 
FS for t'/t=-0.25 at x=0.15 which is shown in Fig.2-c 
(dashed line). At the non-interacting level the change 
between a FS that closes around (tt, tt) and around (0, 0) 
occurs at x=0.22 for the ratio of t'/t used here. In addi- 
tion, in Fig.2-c the experimental points for optimal doped 
YBCO obtained several years ago p2| arc also shown. 



The agreement with our result is only qualitative but 
it has to be considered that the measurements are very 
difficult due to surface effects and thus the experimental 
points have large error bars (not shown) . It is important 
to remark that the similarity found between the interact- 
ing and the non-interacting FS for t'/t=-0.25 does not 
mean that they are identical but the differences will be 
apparent only when larger lattices at lower temperatures 
can be studied. 

Sometimes the experimental results have been inter- 
preted as indicating that there is a FS only along the 
diagonal direction and no FS close to (0, tt) and (71", 0). 
PH The rate of change of n(q) could support this view 
but we found out that it also gets very reduced close to 
(0, 7r) and (tt, 0) even in the non-interacting case when it 
is known that there is a continuous FS. Thus, the present 
results do not allow us to decide one way or the other. We 
also noticed that the interacting FS for t'/t=-0.3 appears 
to be qualitatively different from the non-interacting one. 
In particular at x=0.3 it seems to close around (tt, tt) 
while the non-interacting one closes around (0, 0). Thus, 
if YBCO were to have the same qualitative shape of FS 
than LSCO in the overdoped regime this is another rea- 
son to rule out values of |t'/t| equal or higher than 0.3. 



(0,7t) (7WC) (0, 71) i (Jt,JC) 




(0,71) (7t,7C) 

n — 1 — 1 — 1 — 

\ (c) 



' \ x=0.15- 

- B\ 

(0,0) (7i,0) 

FIG. 2. (a) The numerically calculated Fermi surface 
for U/t=6, (n) = 0.7 on an 8 x 8 lattice for T=0.25t and 
t'/t=-0.25 (closed circles and solid line); the open squares are 
experimental results for LSCO at x=0.3 from Ref. [17]; the 
dashed line is the non-interacting, U/t=0, FS for t'/t=-0.25 
and (n) = 0.7; (b) non-interacting (U/t=0) FS for t'/t=-0.25 
and density 0.9 (dashed line) and experimental data for 
LSCO with x=0.1 from Ref. [17]; (c) non-interacting FS for 
t'/t=-0.25 and density 0.85 (dashed line) together with ex- 
perimental data for YBCO with x sa 0.15 from Ref. [22]. 
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C. Charge Correlations 

The next issue that will be addressed is the origin of 
the incommensurate magnetic fluctuations in the present 
model. In Fig. 3 the charge structure factor N(q) along 
the (0,0) -> (7r,7r) ->• (tt,0) -> (0,0) directions for t'/t 
ranging from (top) to -0.5 (bottom) is shown. In all 
cases there is a broad maximum at Q„ which is due, as 
in the low electron density limit of the Hubbard model, 
to the short range effective repulsion between particles. 

If the incommensurate magnetic fluctuations were due 
to dynamical charge fluctuations, peaks at q = 7r(0, 26) 
and ir(2S, 0) should be observed in N(q) according to 
previous theoretical studies In the present case, 

since 6 — 0.25, the peaks would be expected at (0, n/2) 
and (7r/2, 0). This momentum is indicated with an arrow 
in Fig. 3 and it is clear from the figure that no indications 
of incommensurate charge order is observed. 

Another possible origin of the magnetic incommensu- 
ration in 2D could be simple FS effects. jic|,^,^| There 
are some momenta that map points (or regions) of the 
FS into other points (or regions) also on the FS. These 
are called nesting or pseudonesting vectors and they cor- 
respond to values of momentum where maxima occur in 
the imaginary part of the magnetic susceptibility in the 
2D non-interacting system. We performed a numerical 
calculation of the pseudonesting vector for the interact- 
ing FS. We constructed a histogram in order to identify 
the value of the momentum that mapped most points of 



N(q) 




(0,0) 



(jt.it) 



(Jt,0) 



(0,0) 



FIG. 3. Static charge structure factor for U/t=6, (n) = 0.7 
on an 8 x 8 lattice for T=0.25t and values of t'/t rang- 
ing from (top) to -0.5 (bottom) along the directions 
(0, 0) — (n, 7r) — (-7T, 0) — (0, 0). The arrow indicates the value of 
the momentum where a maximum indicating incommensurate 
short-range order would be expected. 



the FS into other points also belonging to it. For t'/t=- 
0.25 at x=0.3 a maximum in the histogram was obtained 
at q = (n, 0.677r) which is in very good agreement with 
the analytical value for the corresponding non-interacting 
FS, namely q = (tt, 0.7l7r) p3[ . Thus, the maximum in 
5(q) at q = (tt, 0.757r) could be explained by FS effects in 
this case but it may be due to the coarse grid that neces- 
sarily had to be used in our computational studies. In the 
non-interacting case it is expected that the maximum re- 
mains at until x reaches 0.22 (see Sect. III. B). Though 
this behavior seems to be in agreement with previous re- 
sults for the interacting case p(| and in disagreement 
with the experimental data, it is possible that the effect 
of the interaction at smaller dopings will be observed at 
lower temperatures than the ones that can presently be 
reached. 



IV. CONCLUSIONS 

In summary, motivated by recent neutron scatter- 
ing and ARPES experiments, we have studied the U- 
t-t' model numerically on 8 x 8 lattices at temperature 
T=0.25t and with coupling U/t=6. Fixing the density 
at (n) = 0.7 agreement between the position and the 
relative intensity of the incommensurate peaks obtained 
numerically and experimentally for LSCO and YBCO is 
observed for t'/t=-0.25. Larger values of |t'/t| are ruled 
out because in this case the structure factor has a max- 
imum along the diagonal direction rather than at Qg, 
while with less negative values of t'/t a relative max- 
imum is observed at Q„ along the diagonal direction, 
again in disagreement with experimental results. Thus, 
the Hubbard or t-J models without nearest-neighbor elec- 
tron hopping do not reproduce appropriately this exper- 
imental behavior in the cuprates. 

The addition of a diagonal hopping t'/t=-0.25 also pro- 
vides good agreement with experimental measurements 
of the FS for LSCO and YBCO at different densities. A 
FS that closes around (tt, tt) in the underdoped and op- 
timal doped regimes, and around (0, 0) in the overdoped 
case is observed. 

The incommensurate magnetic peaks do not seem re- 
lated to incommensurate charge fluctuations. In the over- 
doped regime FS effects enhanced by the electronic inter- 
actions may be responsible for the observed results. The 
sign problem prevented the exploration of this issue in 
the optimal doped and underdoped regimes. 

The present numerical analysis provides good agree- 
ment between a theoretical model and two unrelated ex- 
periments in the overdoped regime of LSCO. It also shows 
that the consideration of a diagonal hopping in models 
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for the cuprates is crucial in order to reproduce experi- 
mental data. Note that this result is in excellent agree- 
ment with ARPES calculations that have focussed on the 
insulating compound S^CuC^C^. j25| For the cases an- 
alyzed here, the same ratio t'/t can reproduce results for 
both LSCO and YBCO. Although the behavior of the 
incommensurate magnetic peaks in the cuprates appear 
to be similar, there is still no experimental information 
about relative intensities of the peaks at different points 
in the Brillouin zone carried out at the same density. 
These measurements could indicate possible material de- 
pendent properties that could depend on longer range 
electron hopping terms. p5| 
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